Interaction of spin and vibrations in transport through single-molecule magnets 
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We study linear electron transport through a single-molecule magnet (SMM) and the interplay of 
its anisotropic spin with quantized vibrational distortions of the molecule. We show that, despite the 
longitudinal anisotropy barrier and small transverse anisotropy, vibrational fluctuations can induce 
quantum spin-tunneling (QST) and a QST-Kondo effect. The interplay of spin scattering, QST and 
molecular vibrations can strongly enhance the Kondo effect and induce an anomalous magnetic field 
dependence of vibrational Kondo side-bands. 
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Transport measurements on nanometer-sized magnetic 
systems address the fundamental problem of how a few 
magnetic atoms in an anisotropic environment respond to 
an electron current Such an environment is provided, 
for instance, by ligand groups holding such atoms to- 
gether in a single magnetic molecule contacted in a break 
junction [H, A very similar situation arises for trans- 
port through magnetic atoms embedded in a molecular 
network on an insulating surface in an STM setup 
Such systems, for simplicity all referred to as single- 
molecule magnets (SMM), constitute a single, large spin- 
moment with spin-anisotropy. The interplay with quan- 
tum transport provides new possibilities to study and 
control their molecular magnetism. For instance, the 
tunneling allows access to several charge states of the 
SMM which can exhibit enhanced magnetic proper- 
ties (H. When such charge states are only virtually ac- 
cessible, effective spin-spin exchange interaction arises Jj 
and inelastic excitation of the spin moment is possible [2| , 
allowing for time-dependent control [5[. A key result is 
that in either regime the transport depends sensitively 
on the magnetic anisotropy of the SMM, which is char- 
acterized by spin-quadrupole terms in the Hamiltonian. 
A further new aspect is the recently demonstrated me- 
chanical tuning of these anisotropy terms in a transport 
setup [sl. Although the effect of such molecular distor- 
tions on magnetization measurements of SMM crystal 
samples has been addressed [sl, Q , so far their dynamic 
effect on transport through an SMM have not been stud- 
ied. One candidate for sensitively probing such a cou- 
pling of the SMM spin to vibrations is a specific type of 
Kondo effect induced by quantum spin-tunneling (QST). 
This QST through the energy barrier arising from a dom- 
inant uni-axial magnetic anisotropy term relies on the 
presence of weak transverse anisotropy. Combined with 
the exchange tunneling to attached electrodes a QST- 
Kondo resonance specific to SMMs can arise [8]. One 
might expect such QST-assisted Kondo transport to be 
simply suppressed by coupling of the spin to molecu- 



lar vibrations as this tends to increase the anisotropy 
barrier [7]. However, the dynamical effect of vibrational 
fluctuations and the possible competition between longi- 
tudinal and transverse spin- vibration coupling have not 
been studied so far, even though coupling to vibrations in 
the Kondo regime has been considered for spin-isotropic 
molecules (9l-[ll|. 

In this paper we consider the modulation of the mag- 
netic anisotropy of an SMM by a quantized vibrational 
mode distorting an SMM with half-integer spin. Strik- 
ingly, even without static transverse anisotropy^ a QST- 
induced Kondo peak can arise in the linear conductance. 
This Kondo effect is dynamically generated by vibra- 
tional fluctuations which distort the SMM, and thereby 
allow the spin to fluctuate. More generally, a higher 
QST-Kondo temperature may result from spin-vibration 
coupling which is relevant for experimental investigation 
of low temperature transport through SMMs. This en- 
hancement of the interplay of Kondo spin scattering and 
QST by discrete vibrations indicates a possible avenue 
along which transport and quantum magnetism may be 
combined with nano- mechanical effects. 

Model and method. We consider a SMM strongly 
coupled to electronic leads at low temperature in the 
Coulomb blockade regime, where the charge on the SMM 
only changes virtually. We assume that the spin couples 
to a local vibrational mode with frequency ^l. The total 
Hamiltonian reads H = i^sMM + with 



^SMM = 



-DSl + 
-D'St 



JSs. 



(1) 



(2) 



Here Sz is the projection of the molecule's spin on its easy 
axis, which we choose to be the 2:-axis and S± = Sx^iSy. 
We consider here only half-integer values of the spin mag- 
nitude 5, for which there is a Kondo effect at zero mag- 
netic field jsf, and later comment on the integer spin 
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case. Starting from an isolated molecule, the longitudinal 
anisotropy D splits the eigenstates of Sz into the inverted 
parabolic magnetic spectrum, which is sketched in Fig. [1] 
for the representative case of S = 3/2 used throughout 
this work. The zero-field splitting (ZFS) corresponds to 
the energy difference between the ground-state and the 
first magnetic excitation, and equals S = {2S — 1)D for 
^ = (and = D' = 0). The transverse anisotropy 
E breaks the continuous rotational symmetry about the 
easy-axis of the SMM, thereby causing spin-tunneling 
through the barrier. As shown in [7[, vibrational modes 
modulate the magnetic anisotropy and can significantly 
contribute to the observed magnetic splittings. Here 
we additionally consider the dynamical effects of such 
coupling by allowing the dominant anisotropy parame- 
ters to depend linearly on the vibrational mode coor- 
dinate Q = {h ^ b^)/ \f2 through coupling coefficients 
D' and E' . Here the operator h {\)^) relaxes (excites) 
the vibration by one quantum. Thus, when the SMM 
vibrates it lowers its symmetry and QST is enhanced. 
Importantly, this also holds for virtual quantum vibra- 
tions. We note that recently such a linear dependence 
of the D parameter on the pitch angle coordinate [HI 
in the tetra-iron (HI) "propeller" -SMM used in [H has 
been measured. The conduction electron states, repre- 
sented by the operators Ob]^cr->^k(T 0' correspond 
to even combinations of left and right physical electronic 
states and their bandwidth is given by 2VK. Finally, deep 
in the Coulomb blockade regime, the interaction of the 
SMM with the electrodes is given by an isotropic Heisen- 
berg spin-exchange with the conduction band electron 
spin s = \ Y.kk'T.aa'Aa^^^'H'a' ^herc T is the vec- 
tor of Pauh matrices. The coupling J is assumed anti- 
ferromagnetic, which, as pointed out in [Isl, depends on 
the spins of the virtual charge states of the SMM [14 1. 
One might expect that coupling to a molecular vibra- 
tion suppresses the interaction of the SMM with the elec- 
trodes due to Franck- Condon overlap, effectively reduc- 
ing J. However, for an isotropic spin 1/2 it was shown 
that the opposite may happen and that spin-exchange 
processes which change the vibrational quantum number 
are suppressed deep in the Coulomb blockade regime [ol . 
Therefore we assume J to be independent of the vibra- 
tional coordinate in Eq. (j2j). 

The numerical results for the zero-temperature linear 
conductance G(V) = (eVn) (A(ey/2) + A(-eF/2)) as 
a function of the bias voltage V presented below were 
calculated using the Numerical Renormalization Group 
(NRG). Here A{uS) = Y.a ^^(^) is the equilibrium SMM 
spectral function obtained from the T-matrix [15|]. This 
method is exact in the limit of linear response in V or 
of strong asymmetric coupling of the SMM to the elec- 
trodes. For all NRG calculations we used the parameters 
A = 2, TVs = 4000 states [l6| and 11 vibrational states 
were sufficient to obtain results independent of the vibra- 
tion number cutoff. 
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FIG. 1. (Color online). Magneto- mechanical excitation 
spectrum of the SMM with S = 3/2. Magnetic states, ly- 
ing on an inverted parabola due to the static longitudinal 
anisotropy (D), each have vibrational excitations (n) on the 
upright parabolas. Magnetic transitions induced by Kondo 
spin-scattering (J), static anisotropy (E) and longitudinal 
(D') and transverse (E') spin-vibration coupling are indi- 
cated. Magnetic states in the two different Kramer's sub- 
spaces are marked blue and orange. 



Static anisotropy. Even without the vibration {E^ = 
= 0) or transport (J = 0) the SMM eigenstates are 
not spin-eigenstates when E 0. However, for the typi- 
cal case of moderate transverse anisotropy E < D it still 
is convenient to label these mixed states by the dom- 
inant spin-eigenstate (M) in the superposition. Due to 
the two- fold spin-rotational symmetry of Eq. ([T]) the mix- 
ing caused by E is only possible within the two sub- 
spaces spanned by |^|_) (marked blue and orange 
in Fig. [1]). It was shown [sl that upon including ex- 
change spin scattering with conduction band electrons 
(J) the interplay with the QST (generated by E) gives 
rise to a Kondo peak in the linear conductance. Due to 
the presence of the electrodes spin-ffuctuations thus be- 
come significant at low temperature despite the presence 
of the anisotropy barrier of size DS^ opposing SMM spin- 
reversal. A hallmark of this QST-Kondo effect is that it 
is suppressed with decreasing ratio of E/D or increas- 
ing S (because the barrier grows). This QST-Kondo 
effect is clearly distinct from the under-screened high- 
spin Kondo effect which arises for 6* > 1 in the limit 
D = E = D' = E' = Q. = {) where magnetic anisotropy 
is not important [l^. Starting from the latter limit, in- 
troducing the anisotropy barrier, I) > 0, both splits and 
suppresses the high-spin Kondo peak. The remnants of 
the high-spin Kondo peak are located close to the ZFS 
scale 5 = {2S — 1)1), with possible renormalization to 
smaller values for large exchange interaction J Q. 

These ZFS high-spin Kondo side-peaks have recently 
been studied in detail in several experiments [l-Q- The 
QST Kondo peak, on the other hand, located at zero-bias 
in the absence of magnetic field, has to our knowledge 
not been observed experimentally. One possible reason 
for this is that in SMMs typically E/D <1 and the QST 
Kondo temperature Tk is suppressed too much, reducing 
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FIG. 2. (Color online). Logarithmic color plot of the SMM 
spectral function A((jj) normalized to the value Aq — 7r^/4 
achieved for a symmetrically coupled quantum channel. This 
corresponds to the linear conductance at cj = eV and T = 0. 
Parameters S = 3/2, D = 5 ■ lO'^^W , E = O.ID, E' = 0, 
Q = 0.5D, J = 0.2W and D' is varied. The dashed line 
represents the renormalized ZFS from Eq. Q. 



both the height and width of the peak (without splitting 

it). 

Dynamic anisotropy. The anisotropic couphngs of the 
large spin to the vibrational mode and are how- 
ever also of importance especially if the vibrational 
mode frequency energy Q is low. The simplest effect 
of the longitudinal vibrational coupling D' is a pola- 
ronic shift which is different for each magnetic level. For 
E = E' = {) one can shift the vibrational coordinate Q 
(or the operator h) in i^sMM by an S'^-dependent amount, 
resulting in an effective Hamiltonian with eigenvalues 
Em^u = -DW^ - {D''^/n)M^ + ftn, where M = -S,..,S 
and n = 0, 1, .. are the quantum numbers of the spin and 
the polaron, respectively. Each SMM eigenstate is thus 
dressed by vibrational excitations as sketched in Fig. [H 
The effective static energy barrier opposing spin inver- 
sion is changed in shape and increased in height. This 
always results in an effectively enhanced ZFS ^ ^ A, 

A = 5 + {4S{S'-^S + l)-l)^. (3) 

Based on this simple picture one may expect that the 
QST Kondo effect is suppressed by coupling to vibra- 
tions. For fixed E' = this is indeed the case as can be 
seen in Fig. [2j Due to the moderate but non-vanishing 
static E = O.ID a zero-bias QST Kondo peak (red) oc- 
curs which is increasingly suppressed with the coupling 
D' due to an increased barrier. In addition two side-peaks 
are found at the effective ZFS u ^ ±A given by Eq. (j3j) 
as indicated by the dashed lines (and renormalized to 
slightly smaller value due to strong J). These peaks are 
the remnants of the S = 3/2 high-spin Kondo effect. 

In Fig. [3] we now focus on the dynamical effect of the 
spin-vibration coupling by first considering zero trans- 
verse anisotropy {E = 0). Without coupling to the vi- 
brations there is no QST Kondo peak. Switching on spin- 
vibration coupling of only one type, either ^ and 



FIG. 3. (Color online). Effect of the spin- vibration coupling 
on the QST-Kondo peak: SMM spectral functions shown for 
various combinations of zero and non-zero values of the pa- 
rameters. Non-zero values used: D' = 0.04Z), E' = 0.16D 
and E = 0.02D. Remaining parameters as in Fig. O 



E' = (black), or, E' ^ and = (blue), does not 
change this result. Only in the latter case, vibrational 
side-peaks appear in the spectrum at = ±^7 (renormal- 
ized to smaller value due to strong J). Strikingly, when 
both types of couplings are non-zero, a pronounced QST 
Kondo peak appears, even though there is no transverse 
magnetic anisotropy ^ = (red). This vibration-induced 
QST Kondo effect is the central result of this work. We 
now first explain why it requires the presence of both lon- 
gitudinal and transverse couplings referring to processes 
sketched in Fig. [H The Kondo effect is related to fluc- 
tuations between degenerate states of the SMM which 
are in opposite Kramer's subspaces [8]. Since £^ = 0, to 
reach states on the other side of the anisotropy barrier, a 
vibration-induced spin-tunneling of type E' is required, 
which however, involves a virtual vibrational excitation. 
An exchange scattering process (J), which changes the 
Kramer's subspace, cannot change the vibrational num- 
ber at low energy, c.f. This is why the two processes 
E^ and J result only in a QST Kondo side-peak in Fig. [3] 
which is split at Q and suppressed due to the inability 
to reach the vibrational ground-state. Only when a lon- 
gitudinal spin-coupling D' is present as well, the virtual 
vibrational excitation can coherently reach the ground 
state and a full zero-bias QST Kondo anomaly can de- 
velop, as the red curve in Fig.[3]shows. One may say that 
due to the quantum-fluctuations of the vibrational mode 
of the SMM, the magnetic symmetry is broken in virtual 
intermediate states, allowing for Kondo exchange scatter- 
ing, c.f. [lo[. Therefore even for SMMs which have van- 
ishing static E due to symmetry, spin-fluctuations may 
result in pronounced transport features due to the in- 
terplay of exchange scattering and intra-molecular spin- 
vibration coupling. Clearly, the vibrational fluctuations 
can further assist the QST Kondo effect when it already 
is present due to static E: this results in a higher Kondo 
temperature as shown in the green curve of Fig. [3l 

In Fig. m we show the dependence of the QST 
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FIG. 4. (Color online). QST-Kondo temperature Tk in 
units of bandwidth 2W (log color scale) determined from the 
NRG level flow for S = 3/2 and J = 0.2W. (a) Static QST 
Kondo effect: Tk as function of the static anisotropy pa- 
rameters D, E without spin- vibration coupling. Dashed con- 
tour lines indicate that Tk increases with the ratio E/D [3]. 
(b) vibration-induced QST Kondo effect: Tk as function of 
spin- vibration couplings and E' with Q — 0.5D for static 
anisotropy, D = 5-10~^W and ^ = 0. The dashed line marks 
the regime D' > OAD where D' starts to suppress Tk, note 
the constant D here and also the offset: Tk — for either 
E' = or D' = 0. Results for small finite E are similar. 



Kondo temperature Tk on the longitudinal and trans- 
verse static anisotropics and their respective vibration- 
couplings D' and E' . Tk grows as both the transverse 
anisotropics are increased since they both enhance QST. 
Whereas D always suppresses the QST Kondo effect, 
Fig. m^a), its fluctuations first enhance Tk by allow- 
ing QST between the Kramers degenerate ground-states. 
Eventually they will suppress Tk if the vibrational con- 
tribution to the anisotropy barrier S^D''^ /Vt^ protecting 
the SMM against these processes, increases too much. 

Finally we show in Fig. [5] how the experimentally ac- 
cessible evolution of the conductance with a magnetic 
field reveals the different origin of the various peaks. For 
simpHcity, we consider the field to be along the easy 
axis as in some experiments [H, Q and added the Zee- 
man term —SzHz to i^sMM in Eq. ([1]) absorbing the g- 
factor into the magnetic field. As the magnetic field is 
increased, the QST Kondo peak is weakened and splits 
with the anomalous ^-factor, |cjKondo| = 2SHz p^]. This 
clearly indicates the origin of the QST Kondo effect, since 
the ground-state Kramer's doublet M = ±5* is split by 
AM = ±28. Strikingly, the vibrational side-peaks have 
the same, strong field dependence, as they correspond 
to a similar transition offset in energy by Vt\ |cJvib| = 
0.^28 Hz. In contrast to this, the high-spin Kondo peak 
(ZFS) evolves much slower in the magnetic field, inde- 
pendent of the spin magnitude 8: |cjzfs| = A + i^^, 
signaling that it corresponds only to a transition with 
AM = ±1. Comparing with the above formulas for the 
peak evolution, indicated by dashed lines in Fig. O we 
conclude that the QST Kondo and high-spin Kondo ef- 
fects are distinguishable, especially for SMM with large 
spin. Finally, we remark that for an SMM with integer 
spin 8 there is no QST Kondo effect at zero field, but 




^-0.50 



■0.52 



■0.54 



<-0.56 



FIG. 5. (Color online). Magnetic field evolution of the spec- 
tral function (log color scale) for the parameters of the red 
curve of Fig. [S] 5' = 3/2, E = 0.027:>, D' = 0.04L> and 
E' = 0.16D. Dashed lines mark the QST-Kondo (red) and 
vibrational side-peak (yellow) both evolving with an anoma- 
lous g-factor and zero- field split high-spin Kondo peak (blue) . 



instead a zero-bias conductance dip (for < E < D). As 
pointed out in [19] a transverse magnetic field results in 
a QST Kondo effect, where similar spin-vibration effects 
as studied here could occur. 

Conclusion. We have studied the interplay of spin and 
vibration on the conductance through a single-molecule 
magnet. Whereas longitudinal coupling to the vibra- 
tion increases the zero-field splitting, suppressing the 
quantum spin-tunneling Kondo peak, a vibrationally in- 
duced quantum spin-tunneling Kondo effect can occur at 
zero bias if transverse coupling is present as well. The 
transition to virtual vibrational excited states and the 
transverse spin-mixing in these virtual states results in 
a Kondo effect, even in the absence of static transverse 
anisotropy. The interplay with vibrations thus can in- 
crease the quantum spin-tunneling Kondo temperature 
for a given static anisotropy, which may motivate further 
experimental investigation of low temperature transport 
though single- molecule magnets. The measurable mag- 
netic field evolution of the conductance reveals that vi- 
brational side-bands acquire an anomalous ^-factor. We 
acknowledge A. Cornia, J. Kortus and J. Paaske for stim- 
ulating discussions and support from NanoSci-ERA. 
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